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CNJ ' Abstract. We show that nonlinear response of a quantum oscillator displays antiresonant dips and 

resonant peaks with varying frequency of the driving field. The effect is a consequence of special 
symmetry and is related to resonant multiphoton mixing of several pairs of oscillator states at a time. 
We also discuss escape from a metastable state of forced vibrations. Two important examples show 
', that the probability of escape via diffusion over quasienergy is larger than via dynamical tunneling 

■ provided the relaxation rate exceeds both of them. Diffusion dominates even for zero temperature, 
so that escape occurs via quantum rather than thermal activation. The effects can be studied using 
Josephson junctions and Josephson-junction based systems. 

^ : INTRODUCTION 

o 

i-O ' Much progress has been made recently in experimental studies of periodically mod- 

ulated vibrational systems. Examples include optically bistable systems, electrons in 
Penning traps, Josephson junctions, and various nano- and micro-mechanical resonators 
il]-|8]. All these systems display bistability of forced vibrations. Because of thermal or 
^ . externally applied noise, there occurs switching between coexisting stable vibrational 

! states. The measured switching probabilities of noise-induced transitions are in a good 

qh| agreement with the theoretical predictions. 

In the present paper we are interested in the dynamics of a quantum oscillator. The 
generality of the oscillator as a model system and the current interest in quantum 
^ I computing and coherent phenomena lead to two major questions: (i) does a resonantly 

■ driven oscillator display coherent quantum effects that would qualitatively differ from 
those in two-level systems, and (ii) in the presence of relaxation, what is the probability 
of switching between coexisting stable states due to quantum fluctuations? These two 
questions are addressed in the present paper, which is based on the results l^- ITlll . 

A weakly nonlinear oscillator is a multi-level quantum system with nearly equidistant 
energy levels En. Therefore a periodic force of frequency (Of can be nearly resonant for 
many transitions at a time, i.e., hcop can be close to the interlevel distance En+i — En for 
many n. This leads to strong nonlinearity of the response even to comparatively weak 
resonant fields. A well-known quantum effect of the oscillator nonlinearity is the onset of 
Rabi oscillations due to resonant multiphoton transitions. Multiphoton Rabi oscillations 
occur when the spacing between remote energy levels n and m coincides with the energy 
of n — m photons. En — E,n = (n — mjfiCOp fll2l1 . The multiphoton transition amplitude is 
resonantly enhanced, because an m — > n transition occurs via a sequence of virtual field- 
induced transitions k^k-\-l (with m < A; < n — 1), all of which are almost resonant. An 



associated classical effect, in the presence of dissipation, is hysteresis of the amplitude 
of forced vibrations as function of the field amplitude A and cap. 

In this paper (see also [ 10]) we show that multiphoton transitions in the oscillator are 
accompanied by a new effect, antiresonance of the response. When the frequency of 
the driving field adiabatically passes through a resonant value, the vibration amplitude 
displays a sharp minimum or maximum, depending on the initially occupied state. We 
argue that the antiresonance and the multiphoton Rabi oscillations can be observed in 
Josephson junctions. 

If the frequency (Op is close to twice the oscillator frequency, then h(Of is close to 
En+2 — En for many n at a time. This leads to parametric resonance in the oscillator, in 
which it oscillates at frequency (Of /I in response to the driving. Such oscillations are 
intrinsically bistable, because their phase can take on two values that differ by Tt. 

We will be most interested in the semiclassical behavior of the oscillator, which, on 
the one hand, stretches all the way to the classical region, and on the other hand, works 
well for oscillators even deep in the quantum domain. In the semiclassical picture, 
resonant multiphoton transitions correspond to tunneling between Floquet states of 
the oscillator with equal quasienergies. [The quasienergy e gives the change of the 
wave function when time is incremented by the modulation period Xp, \ir{t + 
Xp) — exp(—zeT/7//i) !//(?)]. The occurrence of equal-quasienergy states is related to the 
bistability of forced vibrations of a classical oscillator. 

Tunneling of a driven oscillator is a carefully studiedfis'] example of dynamical 
tunneling [14]. The WKB analysis gives an important insight into the origin of the 
antiresonance, which goes beyond the perturbation theory in the driving field. 

Dynamical tunneling also leads to transitions between coexisting metastable states 
of forced vibrations, which emerge in the presence of dissipation due to coupling to 
a thermal reservoir. In terms of quantum mechanics, dissipation is due to interlevel 
oscillator transitions with energy being transferred to (or from, for nonzero temperature) 
the reservoir. It turns out that dissipation may also lead to transitions between metastable 
states of forced vibrations, even for zero temperature [9]. 

For r = there occur only interlevel transitions where the oscillator energy goes to 
the reservoir (but the energy loss is compensated by the driving field, in the stationary 
regime). However, the quasi-energy may increase or decrease as a result of a coupling-to- 
reservoir induced transition, although with different probabilities. Therefore along with 
drift over quasienergy towards a metastable state, which results from more probable 
transition, there emerges diffusion away from this state as a sequence of less probable 
transitions. The diffusion may lead to activated-like escape. Activation in this case has 
purely quantum nature, and therefore we call it quantum activation. 



THE MODELS 

The Hamiltonian of a nonlinear oscillator with mass M = \ has the form 

H{t) = ^p^ + \(olq^ + \yct^Hp{t). (1) 



We will consider two types of periodic modulation. Hp = Hp'^\ which correspond to 
resonant and parametric driving, 

H^\t) = -qAcos{(OFt), 5(0^'^^ = (Of <Bb, (2) 

H^j''^ = ^q^Fcos{(OFt), Sco'^P^ = ^©F -<»o "C (Oq 

(in what follows we set 7 > 0). 

It is convenient to analyze the dynamics in the rotating wave approximation by 
switching from the fast oscillating operators q,p to slowly varying operators 2,P using 
transformations 

^ = a('")((2cos(oWr + /'sin(oWr), p = -^('"^©^(QsinftjW^-PcosoW?) 
for resonant driving and 

q= a^P\p cos co^Ph-Qsina^Ph), p= -a^P^(O^P\P sin a^Ph + Q cos a^Pk) 
for parametric driving, with 

(oW = ©F, aW = (8©FW'■V37)^/^ 

and 



(0 



= (Of/2, a(^) = (2F/37)V2. 



The variables Q,P arc the appropriately scaled coordinate and momentum. The com- 
mutation relation for them has a simple form 



[P,Q]^-a, A = '(a^'"'^)) ^ (3) 

The parameter A plays the role of the effective Planck constant. We note that it is 

proportional to the oscillator nonlinearity 7 scaled by the comparatively small detuning 
of the field frequency, in the case of nearly resonant driving, or the comparatively small 
field amplitude, in the case of parametric driving. 
The dynamics of Q, Pin the two cases is described by effective Hamiltonians 

h(p^ = {F /4) (a^P^^ g^Pl (4) 

Their eigenvalues are equal to the quasienergies e„ of the oscillator. The functions g^^'P^ 
have the forms 

g^'\m^\{Q^+P^-lf-fi''^Q. j8 = (5) 

gip) = ^{Q^ + pY + ^-{l-li)P^ -1(1+^)2', M = 2cOFdco^P^/F. (6) 



They are shown in Fig. \T\ Each of them depends on one parameter. In the region of 
bistability of period one vibrations, < /3 < 4/27, the function g^''^ has a shape of a 
tilted Mexican hat, with a maximum at the top of the central dome and a minimum at the 
lowest point of the rim. For a parametrically excited oscillator in the region — 1 < /i < 1 
the function g(^) has two symmetrical minima. These extrema of correspond to 
metastable states of the driven oscillator in the presence of weak dissipation. The saddle 
points of g('''^^) correspond to unstable states of forced vibrations. 




FIGURE 1. The scaled classical quasienergy of the oscillator for resonant and parametric driving (left 
and right panels, respectively) as function of the slowly varying coordinate and momentum Q,P. The plot 
of g'^''^ refers to the reduced field intensity p = dyA^ / 32 0)j {5(0^''^)^ = 1/27. The plot of g^P^ refers to 
H^2cof5co'^p^/F = -0.1. 



MULTIPHOTON ANTIRESONANCE 

We will start with the studies of the coherent response of the oscillator to a nearly 
resonant field. When the field amplitude A 0, the eigenstates \n) of the Hamiltonian 
coincide with the Fock states of the oscillator, and the quasienergies are 

£„ = -5(0^/1 + Vn(n+1), V = ^. (7) 
2 4(Uq 

We keep only the lowest-order term in V, which corresponds to the weak nonlinearity 
approximation. In this approximation the energy of an Mh oscillator state for A = is 
E/^ = hcooN + VN{N + l)/2. The A/^-photon resonance A/^^CQ^- = Ej^ — Eq occurs, in terms 
of 5(0^, for 

5(oW = 5£oi,'"^=V(A^+l)/2. 

For the corresponding field frequency (Of the quasienergies Cq and are equal. 

The field leads to mixing of the wave functions of resonating states and to level 
anticrossing. This anticrossing is clearly seen in the upper right panel of Fig. |2| The 
minimal splitting of the levels Cq and is given by the multiphoton Rabi frequency Q.r. 
For weak field it can be obtained by perturbation theory 1 12]. To the lowest order in the 



field amplitude A in the limit of large A'^ 

aR = ViA/ANfN^/\2n)-^/\ A^ = {2h(Oo)^/^\V\N^/^ £xp{-3/2)/2. (8) 

The Rabi frequency depends on exponentially, Q.r oc . This dependence works well 
numerically in the whole region A < Ayy 1 10]. 




FIGURE 2. Quasienergy levels e„ and reduced susceptibilities in the resonating states of a driven 
nonlinear oscillator The plots refer to a 5-photon resonance, N = 5. The labels 1 and 2 correspond to 
the ground (n — 0) and 5th (n ~ 5) oscillator states in the absence of driving, for 5(0 = Sco^'^ < SOj'^'. 

Left panel: the Umit of zero driving amplitude, A ^ 0. The quasienergy levels cross for 5(o)^ /V = 3. 
The susceptibilities in the resonating states also cross at that same frequency. Right panel: level repulsion 
and antiresonance of susceptibilities for comparatively weak modulation. The data refer to A/An = 1/2. 
It is seen from the figure that the susceptibility antiresonance is much more pronounced than the level 
repulsion. The scaling factor for the amplitude a„ of oscillator vibrations in the states « = 0,5 is / = 
{Mcoo)-^/^A. 

Coherent response of the oscillator to the driving field is characterized by the expecta- 
tion value of the coordinate q. If the oscillator is in an eigenstate \n) of the Hamiltonian 
H^''\ this value has the form 

qn = {h/2cOoy/^a„e-''^^' +C.C., q„ = {n\q\n). (9) 

To first order in the field, the reduced amplitude of forced vibrations in an nth state an 

is 

a„ = -/5(D(''V((5(oW-Vn)[5(oW-V(n + l)]), f = {Shcoo)-^/^A. (10) 

Remarkably, for 5(0^'^^ = Scoj;^^ the vibration amplitudes in the resonating states coincide 
with each other, ayv_„ = a„ for < n < N/2, see the left lower panel in Fig. El 

Field-induced multiphoton mixing leads not only to splitting of the quasienergy levels, 
but also to repulsion of the vibration amplitudes. It can be calculated by diagonalizing 



the Hamiltonian H^^^ and is shown in the right lower panel of Fig. |21as a function of 
frequency detuning 5(0 = 5(0^'^\ One of the involved resonating states is the ground 
state of the oscillator n = in the limit A 0. The quantities plotted in Figs. Etb) 
are susceptibilities, they are proportional to the ratio of the vibration amplitude to the 
modulation amplitude a„/A. 

The antiresonant splitting of the expectation values of the vibration amplitudes is by 
far the most interesting feature of Fig. El It occurs at the adiabatic passage of 5(0^'^^ 
through resonance, where the system switches between the ground and excited states. 
In particular, the amplitude displays an antiresonant dip if the oscillator is mostly in the 

ground state for (5(o('') - 5(ol^^)/V < 1 or in the state for (5 - 5(oj^^)/V > 1. The 
magnitude and sharpness of the dip are determined by Q.r/V and depend very strongly 
on the field and A^. With decreasing Q.r/V the dip (and peak) start looking like cusps 
located at resonant frequency. We note that, in contrast to the case of energy levels, there 
is no reason for repulsion (anticrossing) of susceptibilities. In fact, as seen from Fig.|21 
the susceptibilities do cross, although away from the resonant frequency. The effect of 
antiresonance is due purely to specific quantum interference [ 10]. 

The dip in the oscillator response has no analogue in two-level systems. There, for 
nearly resonant driving, the coherent response in the two adiabatic states differs only in 
sign. It displays a peak when the radiation frequency adiabatically passes through the 
transition frequency. 

The WKB picture of the antiresonance 

In the WKB approximation, Rabi oscillations correspond to tunneling between the 
states with nearly equal quasienergies. Such semiclassical states can be found from the 
Hamiltonian H^''^ ^ by using the Bohr-Sommerfeld quantization condition applied to 
the mechanical action §PdQ for trajectories g{Q,P) = const, with h replaced by A, 
Eq. 01). It is seen from the left panel of Fig.[Tlthat, in a certain range of g, there are two 
types of trajectories with the same g, those on the internal dome and those on the external 
part of the Mexican hat g{Q,P). If, as a result of the Bohr-Sommerfeld quantization, the 
quantized values of g on the two parts of the surface g{Q^P) coincide with each other, 
then there may occur resonant tunneling between the corresponding quantum states. The 
resulting tunneling splitting llisll is the Rabi frequency. 

Interestingly, one can show that the average value of the coordinate 

Q{g) = ^'\g) / Q{t)dt 
Jo 

[t{g) is the period of oscillations for a given g] is the same for the internal and external 
trajectories with the same g. This corresponds to the susceptibilities of the resonating 
quasienergy states being the same, in the neglect of tunneling-induced mixing of the 
states. We emphasize that the fact that the susceptibilities are equal is not a result of the 
perturbation theory in the field amplitude, as in the case of Eq. (flOl) . they are equal in 
all orders of the perturbation theory in A as long as tunneling is disregarded. Tunneling- 
induced state mixing leads to the antiresonance of the response lHoi] . 



ESCAPE OF A DRIVEN SYSTEM: TUNNELING OR QUANTUM 

ACTIVATION? 



We will now briefly outline the new effects and unanswered questions that emerge when 
dissipation is taken into account 10, [HI]. We will assume that, even though dissipation 
is weak, the dissipation rate exceeds the tunneling rate. The problem of fairly general 
interest that will be addressed is switching between metastable states of forced vibrations 
of a quantum oscillator. We will consider the most interesting situation where there are 
many quasienergy states between the extrema of the quasienergy surface. In the case 
of escape of a particle from a potential well it corresponds to a well with many energy 
levels. 

In systems in thermal equilibrium, the rate of tunneling decay of a metastable state for 
low temperatures is given by the probability of a tunneling transition from the ground 
state in a metastable potential well. In the case of a resonantly driven oscillator this 
corresponds to dynamical tunneling from the top of the dome of the quasienergy surface 
g^''^ to the state on the external orbit with the same quasienergy, see the left panel 
of Fig. [U The tunneling is shown schematically in the central panel of Fig. 13 For a 
parametrically driven oscillator the corresponding tunneling occurs between the minima 
of the surface g^P\ as shown in the right panel of Fig.|3l 

For higher temperatures, again in the case of equilibrium systems, one has to take 
into account tunneling from excited intrawell states. Escape may occur also via ther- 
mal activation over the potential barrier. One of these escape mechanisms dominates, 
depending on temperature [15]. In the case of a driven oscillator tunneling from excited 
states corresponds to tunneling with quasienergies that differ from those at the extrema 
of siQ^P)- In addition, there is a probability of activation over the quasienergy barrier. 
However, since the distribution over quasienergy is not of the Boltzmann form, it is not 
clear which of the escape mechanisms dominates at a given temperature of the bath. 

For weak dissipation the distribution can be described by the balance equation for the 
occupations p„ of quasienergy (Floquet) states \n), 

Pn = - Wnmpn + W,nnpm- ( 1 1 ) 

The transitions probabilities W„m can be calculated as matrix elements of the operator 
that describes relaxation of the oscillator. The wave functions \n) can be found from the 
Bohr-Sommerfeld approximation disregarding tunneling. A standard WKB calculation 
allows one to express Wnm in terms of the Fourier components of the coordinate and 
momentum Qm-n{gn),Pm-n{gn) for a given quasienergy gn [the functions Q{t),P{t) 
are periodic functions of time for a given g, with period T{g)]. In the simple case of 
linear friction, which corresponds to relaxation transitions between nearest energy (not 
quasienergy) levels of the oscillator shown in the left panel of Fig. |3l Wnm are simply 
quadratic in Qm-n{gn)-,Pm-n{gn)- They exponentially decay with \n — m\. 

The probabilities W„,„ are organized so that it is more likely for a system to make 
a transition toward the value of g in the metastable state rather than away from it. 
This is why the state is metastable. However, in contrast to systems in equilibrium, the 
probabilities Wnm do not satisfy the condition Wnm = ^mn^'^vliSn — gm) /kT]. Even for 
r there is a nonzero probability to make a transition in the direction opposite to the 



metastable state. This is a consequence of the fact that the Floquet states \n) are linear 
combinations of the Fock states of the oscillator. Therefore, even where all transitions 
between the Fock states go in one direction in energy, as in the left panel of Fig. El 
transitions between the Floquet states go in different directions in quasienergy, although 
with different probabilities. 
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FIGURE 3. Left panel: transitions between the energy levels of the oscillator £„ due to coupling to a 
thermal reservoir. For T ~ Q there occur only transitions to lower levels with emission of excitations of 
the reservoir. In the simplest case relaxation is due to transitions between neighboring oscillator levels. 
Central panel: the cross-section of the surface of the scaled quasienergy ' shown in Fig.^by the plane 
P — 0. Quasienergy levels for orbits that lie on the dome in Fig.^are shown schematically. Relaxation to 
the metastable state at the top of the dome corresponds to transitions to higher quasienergy, whereas less 
probable transitions to lower quasienergy lead to an activation-type escape from the metastable state. Right 
panel: the cross-section of the quasienergy surface of the parametrically driven oscillator g^P\ Fig.^ by 
the plane P — Q, with schematically shown quasienergy levels. Relaxation drives the oscillator down to the 
lowest quasienergy levels at the bottom of the wells of g^'^\Q,P), but there is a possibility of activation 
even for T = 0. 

Transitions in the "wrong" direction lead to diffusion over quasienergy away from 
the metastable state. Their immediate consequence is a finite width of the stationary 
distribution over quasienergy even for T — > 0. Another closely related consequence is 
a nonzero probability to reach quasienergy of the saddle of the functions starting 
from a metastable state. The saddle of g^'''P^ is similar to the top of a potential barrier 
for a particle in a metastable potential. The logarithm of the probability of reaching the 
saddle gives the activation energy of escape as a result of diffusion over quasienergy. We 
call this process quantum activation, since it occurs even for T = 0. 

For resonantly and parametrically driven oscillators the transition rates W„m can be 
calculated explicitly using the fact that the classical trajectories in (2,/') -variables can 
be expressed in terms of the Jacobi elliptic functions [9, 11]. We have compared the 
activation energies of escape with the tunneling exponents. When the activation energy is 
smaller than the tunneling exponent in the absolute value, escape occurs via quasienergy 
diffusion. We found that, both for a resonantly and a parametrically excited oscillator, 
escape from metastable states occurs via activation, not tunneling. This holds true for all 
parameter values where a nonlinear oscillator has coexisting stable states. We found an 
unusual behavior of the distribution for an underdamped parametrically driven oscillator 
for r = 0. For some quasienergy the distribution displays a sharp decrease, and at the 
same time the tunneling rate goes to zero. 



The physical origin of the fact that escape occurs via activation, not dynamical tun- 
neling, remains not fully understood, the existing arguments are formal [9]. Apparently, 
there must be a crossover from activation to tunneling when the system goes to equi- 
librium, but the models that we have discussed are strongly nonequilibrium, the very 
presence of coexisting metastable states is due to periodic driving. 

In conclusion, we have shown that a simple system, a driven nonlinear oscillator, 
displays unusual quantum coherent phenomena and unusual switching behavior. The 
considered effects have no analogue in two-level systems and are qualitatively different 
from what has been known about switching in thermal equilibrium systems. They are not 
only of fundamental interest, but are also important for many applications, in particular 
in sensing and quantum measurements. 
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